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H \ strong coupling, large-N expansion. In this limit we explicitly perform the func- 

tional integral over fundamental degrees of freedom and obtain the effective chiral 
Lagrangian. We perform this calculation on the Body Centered Hypercubical lat- 
tice which preserves Euclidean invariance to order p^. We further discuss how the 
coefficients could be obtained numerically, out of the strong-coupling domain. 
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1 Introduction 

Effective chiral Lagrangians, originally introduced by Weinberg |[1[ , provide a very power- 
ful method to describe the low energy phenomenology of QCD, exploiting the consequences of 
chiral symmetry and current algebra 0. The idea underlying effective Lagrangians is simple 
and elegant. There are many phenomena which, although the fundamental theory involves 
several different scales of length and energy, only lead to the excitation of a subset of all the 
degrees of freedom. We will call such degrees of freedom the effective degrees of freedom. 
Under these circumstances it is conceptually possible to integrate out of the theory all of the 
other degrees of freedom and describe the phenomena in terms of an alternative Lagrangian, 
namely the effective Lagrangian, more amenable to calculations. Thus in low energy QCD, 
below the threshold for production of resonances, the effective degrees of freedom are those of 
the pseudoscalar fields, i.e. the would be massless Goldstone bosons of spontaneously broken 
chiral symmetry, and pion and kaon interactions can be described in terms of an effective 
chiral Lagrangian. Another example is found in the physics of heavy quarks, where effective 
non-relativistic Lagrangians can be used to describe a large class of processes which do not 
lead to the excitation of the relativistic degrees of freedom of the heavy quarks. 

The predictive power of effective Lagrangians in general derives from the existence of some 
underlying exact or approximate symmetry, such as chiral symmetry in QCD, which restricts 
substantially the type of couplings that can appear in the Lagrangian itself. Moreover, 
these couplings can be arranged in an expansion into powers of some parameter related to 
the suppression of the non-effective degrees of freedom (e.g. the magnitude of the external 
momenta in QCD) and at any order in the expansion there is only a finite number of possible 
coupling terms. In typical phenomenological applications, then, the corresponding coupling 
constants are fixed by comparing the results that follow from the effective theory to a subset 
of the experimental data. This determines the effective Lagrangian, which can subsequently 
be used to calculate a variety of other experimental observables 0, ^ . 

Although this way of proceeding provides a very powerful computational tool, capable of 
producing a large number of theoretical predictions starting from general symmetry consider- 
ations and a limited amount of experimental information, it would be even more valuable to 
be able to calculate the coefficients of the effective Lagrangian directly from first principles. 
In the case of low energy QCD, this goal appears within reach of advanced lattice numerical 
simulations and in this article we wish to present a framework by which it can be achieved. 

In some sense, once one has the possibility of calculating observables by means of numer- 
ical simulation methods, deriving the coefficients of the effective Lagrangian might appear 
straightforward: one would just need to compare a sufficient number of predictions from the 
effective Lagrangian to the results of numerical calculations, very much as one does compar- 
ing them to experimental data in the phenomenological derivation. But then, from this point 
of view, why resort to an effective Lagrangian at all, if the observables can be calculated di- 
rectly? The point is that numerical simulations provide only an approximate computational 
tool, where the accuracy of the calculations is severely limited by statistical errors and the 
discretization procedure. The precision by which one can determine quantities such as the 
coefficients of an effective chiral Lagrangian, or even the feasibility of such a calculation in 
absolute, depends dramatically on the choice of a sound computational technique. 



The basic idea underlying the method we wish to propose is that the separation into 
effective and non-effective degrees of freedom of the continuum theory will be mirrored by a 
similar distinction in the theory discretized over a lattice. It will then be possible to define 
a lattice effective Lagrangian to be used as an intermediate step in the derivation of the 
coefficients of the continuum chiral Lagrangian from the fundamental theory discretized over 
the lattice. The advantage of this two step procedure is that the lattice effective Lagrangian 
will contain explicitly the collective fields responsible for the long distance behavior of the 
fundamental lattice theory. The calculation of its coefficients, which for realistic applications 
will have to be made through a numerical integration of the non-effective degrees of freedom, 
will then only involve shorter scales of length, and should be feasible on a lattice of moderate 
size. Being able to reduce the size of the lattice, and therefore the total number of variables 
in a numerical simulation, produces major advantages for the accuracy and efficiency of the 
calculation. 

In this article we will illustrate these ideas by deriving the lattice effective chiral La- 
grangian in the strong coupling approximation and using it to calculate, always for strong 
coupling, the coefficients of the p"^ terms of the continuum Lagrangian. This will be done in 
Section 3, after a quick recapitulation, presented in Section 2, of the form of the continuum 
chiral Lagrangian and of some features of the lattice discretization of the QCD fundamen- 
tal Lagrangian that will be relevant for the rest of the paper. Section 3 owes a lot to the 
truly pioneering works of Kluberg-Stern, Morel, Napoly and Peterson ^ and of Kawamoto 
and Smit [|] where the methods for integrating over quark and gauge degrees of freedom in 
the strong coupling limit were clearly laid out and an effective strong coupling Lagrangian, 
expressed in terms of meson fields, was first derived. We add to their work, however, the 
consideration of a different lattice, namely the Body Centered Hypercubical (BCH) or F4 
lattice, whose greater symmetry implies that invariance under lattice transformations carries 
over to Euclidean invariance up to the terms of order p"^ in the continuum. This is clearly 
crucial for a strong coupling derivation of the p^ terms in the effective chiral Lagrangian. We 
were inspired to use this lattice by a talk presented by Boghosian at the Amsterdam Sympo- 
sium on Lattice Gauge Theories. This talk ^ focused on the application of the BCH lattice 
in the theory of cellular automata and lattice gases, however the p^ invariance had been 
previously emphasized and exploited by Neuberger |p in a study of the Higgs phenomenon. 
In the context of lattice field theories, this lattice was also considered in references [^ - [pT| . 

Finally, in Section 4, we will discuss how one could derive the parameters of the chiral 
lattice Lagrangian for more realistic values of the QCD coupling constant, closer to the 
continuum limit, by numerical simulations coupled to real space renormalization techniques. 

2 Chiral Perturbation Theory and Lattice QCD 

For the convenience of the reader and in order to establish a point of reference we review 
the basic features of chiral perturbation theory^ and of the lattice formulation of QCD. 

Following Gasser and Leutwyler [Q and ^, we introduce the generating functional for 
the connected Green's functions of QCD in the presence of scalar, pseudoscalar, vector and 

^Several reviews are given in references [12| - nm. 



axial vector sources^ 

^w[3,P,vM =< I Te'I'^^''^ I >, dl 



^ = ^QCD - qix)[s{x) - i-f5p{x)]q{x) + g(x)7^[t;^(x) + 75a^(x)]g(x). (2) 

In this expression, q, q represent the hghtest quark flavor triplet and CPqcd is what remains 
of the Lagrangian of QCD when the masses of three lightest quarks are set to zero 
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^IcD = -—,G,,{x)G^^{x) + q{x)^^[id, + G,{xMx). (3) 



s(x),p(x), f^(x) and a^(x) are Hermitian matrices in flavor space. The quark mass matrix 



rris ) 



(4) 



is contained in the scalar field s{x). In this paper we will not consider the ^[/(S) singlet 
vector and axial vector current: 

Ti\v^{x)\ = Tr[a^(x)] = 0. (5) 

It is convenient to define the right and left handed currents: 

r^(x) = v^{x) + a^(x), l^{x) = v^{x) - a^(x). (6) 

The Lagrangian (|^) exhibits a local fiavor chiral symmetry SUl{^) x SUh^S): 

q{x) -> L{x) ^ ^ q{x) + R{x) — ^^qix), 

ffiix) — > R{x)r^{x)R{x)'^ + iR{x)d^R{x)\ 
l^{x) -^ L{x)l^{x)L{xy + iL(x)dfj,L(xy , 
s{x) + ip{x) -^ R{x)[s{x) + ip{x)]L(xy , (7) 

with matrices L{x),R{x) G SU{3). 

The goal of chiral perturbation theory is to produce an expansion of the generating 
functional W in powers of the external momenta and quark masses. Because of the presence 
of pseudo Goldstone bosons one can not simply perform a Taylor series expansion of W, but 
one must introduce suitable fields to account for their degrees of freedom. These fields are 
collected in a unitary 3x3 matrix: 

U{x)=exp[^^^^^], U{x)U{x)^ = l. (8) 

A", a = 1, ..., 8 are Hermitian generators of SU{3), with the normalization 

Tr(A"A^) = 2(5'^^ (9) 

^Although it is in principle possible to include tensor sources, they are usually neglected. 



7r"(x) are real pseudoscalar fields and -F° is the pion decay constant in the chiral limit: 

A" 

< I g(a;)7^75Yg(x) | vr* > = if^p^6''\ 

/^ = Fo(l + 0{mq)) ^ 93.3MeV. (10) 

To the lowest order O(p^), one must only evaluate tree level diagrams with a classical La- 
grangian involving terms with two derivatives and one power of scalar and pseudoscalar 
sources: 

£2 = -fTT{[D^U{x)][D^'U{x)]^ + x{x)U{xy + x(x)^f/(x)}. (11) 

In this expression, 

D^U{x) = d^U{x) - ir^{x)U{x) + iU{x)l^{x) (12) 

is a flavor covariant derivative and matrix x(x) collects the scalar and pseudoscalar fields: 

X{x) = 2Bq[s{x) + ip{x)]. (13) 

The constant Bq is related to the quark condensate in the chiral limit: 

< I Mu I >o=< I rfd I >o=< I ss I >o= -F^Bo{l + 0{mq)). (14) 

In the expansion in powers of pion momenta, vector and axial vector sources count as contri- 
butions of order p, whereas scalar and pseudoscalar sources count as contributions of order 

To the next to leading order O(p^), one encounters two kinds of contributions: 

• one loop diagrams with £2 vertices; 

• tree diagrams with £4 vertices. 

£4 is the most general Lagrangian which contains four powers of momenta consistent with 
the local chiral flavor symmetry (^). Therefore, it can be formed only from the covariant 



derivatives (p!2| ) and the fleld strengths for right and left-handed flelds: 



F^" = 9^r -9^P -i[P,r]. (15) 

£4 = L,{TT[{D^U){D^U)^y + L2TT[{D^U)\DM)m{D^U)\D''U)] 

+ LsTi[iD^U)\D^U)iDM)HD''U)] + L,Tr[iD^U)\D^U)]Tiix^U + xU^) 

+ L,Tt[{D,U)\D^U){x^U + xU^)]+L,[Tt{x^U + xU^)? 

+ LrlTiix^U - xU^)? + LsTiix^Ux'U + xU^U^) 

- tLgTi[Fll{D''U)iD''Uy + FI:^{D>'U)\D''U)] + Liq1:i{U^ F^^U F^^"") 

+ i7iTr(F^^,F^^^ + F^^X'^'^) + i72Tr(x^x). (16) 

As we see, at leading order, chiral symmetry restricts the number of terms in the effective 
chiral Lagrangian to only two, so that one needs to know only two constants, Fq and i?o, in 



order to describe all low energy phenomena in QCD to order p^. At the next to leading order 
one needs ten more constants: Li - Liq. (The terms with coefficients Hi and if 2 are of no 
physical significance and we will disregard them from now on.) 

In two seminal papers by Gasser and Leutwyler ^, ^ these ten coefficients have been 
calculated to reasonable accuracy by comparison with experimental data. It would be of 
course important to be able to determine all the coefficients theoretically. This requires 
a non-perturbative calculation and lattice numerical simulations represent today the only 
method by which such a calculation could be carried out from first principles. Otherwise one 
can try to calculate the coefficients from models which are similar to, but easier to solve than 
QCD. For an account of various attempts at a theoretical determination of the coefficients 
see, for example, []T5|, |T6|, |I^ and references therein. In recent years, a rather exhaustive 
calculation has appeared in references [l^, |l^ . 

Our goal is to provide a framework by which the coefficients of the terms of order p^ in 
the chiral Lagrangian can be efficiently calculated by means of a lattice numerical simulation. 
The idea is to proceed through the intermediary of an effective chiral Lagrangian defined over 
the lattice. The problem then splits into two separate calculations: 

• A numerical simulation is used to calculate the coefficients of the lattice chiral La- 
grangian through a matching of lattice observables. As we will argue in Sect. 4, it 
should be possible to carry out this step on a reasonably small lattice and, therefore, 
with high statistical accuracy. 

• The coefficients of the continuum chiral Lagrangian will be obtained from those of 
the lattice effective Lagrangian by an expansion of the observables for small lattice 
momenta and by carrying out the appropriate renormalizations. This step can be 
performed analytically, by a perturbative expansion of the lattice chiral Lagrangian, if 
the coefficients, as determined in the previous step, are such as to validate the expansion. 
Otherwise it would have have to be done numerically. 

We believe that one can carry out these two steps with better accuracy than a direct 
determination of the coefficients of the continuum Lagrangian through a comparison with 
observables calculated over the lattice. 

While in Sect. 4 we plan to discuss in some detail the procedures that could be followed to 
perform such calculations in the physically interesting domain where continuum scaling sets 
in, the main result which we wish to present in this paper is an illustration of the method in 
the strong coupling approximation. Precisely, we wish to establish the possibility of defining 
a lattice chiral Lagrangian by deriving it explicitly at leading order in the strong coupling 
expansion in /3 = Q/qq- An expansion for low lattice momenta will then allow us to identify 
the values of the coefficients of the continuum Lagrangian. 

In the balance of this section we review a few basic features of the lattice regularization 
of QCD that will be useful for our subsequent derivation of the lattice chiral Lagrangian and 
will also serve to establish our notation. 

In the lattice regularization of QCD, Euclidean space-time is generally discretized by the 
introduction of a regular lattice of points and oriented links joining neighboring points. In 
the vast majority of applications this lattice is a hypercubical lattice, however, for reasons 



that will become clear in Sect. 3, we will need to work with a different lattice to carry out 
the strong coupling expansion. Therefore we will leave the detailed geometry of the lattice 
unspecified for now, but we will assume that the lattice is regular and uniform, with lattice 
spacing a. 

We will denote by x the coordinates of the lattice points and by v the fundamen- 
tal displacement vectors in the lattice. Thus, for clarification, with a hypercubical lat- 
tice we would have x = {iia,i2a,i^a,iia) with integer ii,i2,i3,i4: and v would range over 
(a, 0, 0, 0), (0, a, 0, 0), (0, 0, a, 0), (0, 0, 0, a). 

The fundamental variables are the gluon fields and the quark fields. The gluon fields are 
finite elements of color SU{3), which are associated to the oriented links of the lattice. It is 
rather standard to denote such variables by Ux,v, but eventually, following another standard 
notation, we will want to reserve the U symbol for the chiral fields. We will therefore denote 
the lattice gauge fields by G^^v 

The gauge part of the action is constructed from the plaquette variables, i.e. the transport 
factors around the elementary closed contours of the lattice. The plaquettes are defined by 
two displacement vectors, and so we will use the notation Gx,viV2 to denote a plaquette 
variable. (For the hypercubical lattice Gx,viV2 = Gx,viGx+vi,v2Gx+v2,viGl.^,J. The gauge part 
of the action is then given by 

Sg = l3K, J2 Tr(/-G,,,,,J, (17) 

where we are summing over all plaquettes with both orientations, Ki, as well as K2 in the 
equation below, are lattice dependent normalization factors, and the coupling parameter (3 
is related to the bare coupling constant by /5 = Q/qq- 

The quark and antiquark fields J^ ipx,'ipx are defined over the sites of the lattice. The 
interaction term between quarks and gluons is given by: 

S^ = a^K2 Y^ Y^{i)xi}Gx,vi^x+v - i^x+vi'Gl^^ipx) (18) 

X V 

We will also include scalar and pseudo-scalar source terms through a matrix jx'- 

Sj = «^ H i^xjx'ijJx = a'^Yl ^^(^^ ~ ilbPx)tpx- (19) 



To summarize, equations (p^pJ| ) define the QCD action on any regular uniform lattice. 



3 Strong coupling expansion on the BCH Lattice 

In order to derive a lattice chiral Lagrangian one must integrate out the high energy 
degrees of freedom, i.e. the quantum fluctuations of the gauge and quark fields, re-expressing 
the generating functional in terms of pseudoscalar fields conjugate to the external sources. 
Thus one effectively performs a bosonization of the original theory. For general values of the 



"^Ki = K2 = 1/2 for the hypercubical lattice. 

^ Color, spin and flavor indices will be left implicit where it doesn't affect the clarity of the formula. 



bare lattice coupling constant this can only be done numerically, but in the strong coupling 
and in the large A^ (number of colors) limits the bosonization can be done analytically. The 
techniques for deriving an effective Lagrangian in the strong coupling limit were laid down 
in papers by Kluberg-Stern, Morel, Napoly and Petersson [Q and Kawamoto and Smit [^]. 
We add, however, to the work of these authors an ingredient that plays a crucial role for the 
possibility of obtaining the terms of order p^ in the chiral Lagrangian. As in most lattice 
calculations, P, H deal with a hypercubical (HC) lattice. While such a lattice obviously does 
not have the symmetry of the continuum, all tensors of rank up to two symmetric under the 
group of lattice transformations are also symmetric under the full group of 4-dimensional 
rotations of Euclidean space-time. However, this property does not carry through to tensors 
of higher order, so that, although one could expand the chiral Lagrangian derived on a HC 
lattice to terms of 4^^ order in lattice momentum, one would obtain terms that cannot be 
identified with continuum counterparts. This precludes the possibility of recognizing the 
structure of the continuum chiral Lagrangian to 0{p'^) on the HC lattice. 

In order to pursue our calculation, we will study QCD on a different lattice, namely 
the Body Centered Hypercubical (BCH) or F4 lattice, which, as already mentioned in the 
introduction, has a remarkable property that not only tensors of rank two but also those of 
rank four have the continuum symmetry. 

The BCH lattice can be formed in two equivalent ways. One is to visualize it as the usual 
HC lattice from which all the sites with I]^x^ = odd (or even) have been removed. An 
equivalent way is to take all the sites of the HC lattice together with centers of its elementary 
cells (hence the name BCH) to represent the sites of the new lattice. The lattice formed in 
such a way has the largest symmetry group of all four dimensional regular lattices. Whereas 
the usual HC lattice has 384 element symmetry group, the symmetry group of the BCH 
lattice has 1152 elements. A consequence of this is the required p^ invariance. 

On the BCH lattice every site has 24 nearest neighbors forming triangular plaquettes. In 
terms of unit vectors of the hypercubical lattice Cj, the links in the positive direction f "• are 
given by the following twelve linear combinations: 



O7 I Cx O /J 



<. = ^^^^^,{l<.<j<4,a = ±l}. (20) 



The use of the BCH lattice entails another advantage in the strong coupling limit. It is 
of rather technical nature and so we will mention it only briefly. The definition of lattice 
fermions encounters some notorious problems. A straightforward discretization of the Dirac 
operator leads to the introduction of additional poles in the propagators of the fermions at the 
corners of the Brillouin zone (species doubling). Although these extra fermionic modes would 
in general be coupled through the gluonic quantum fluctuations, on a hypercubical lattice 
there is a fourfold degeneracy of totally decoupled fermionic degrees of freedom. This is akin 
to starting with one flavor of fermions, but discovering that the theory actually contains 
four decoupled flavors. As a consequence, the spontaneous breaking of chiral symmetry gives 
origin to 16 pseudo Goldstone bosons, even though one started with a single flavor. On a 
BCH lattice, because of the higher coordination number of its vertices, such a separation of 
the fermionic fields into four decoupled components (which is not demanded by any general 
theorem on lattice fermions) does not take place and one maintains the relationship between 
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the number of explicit flavor indices and the number of pseudo Goldstone bosons proper of 
the continuum theory.^ 

In order to study interactions with external vector and axial vector fields, we will introduce 
traceless Hermitian sources on the positive link v starting from a lattice point x: 

^).,„ = <,A^ S,,„ = <,,A^ (21) 

We couple these sources to quark and gluon fields in a nonlinear way through an auxiliary 
fieldy({),fi): 

y{v, a):,., = exp{- — [v^^^ + v^+^^^ + 75(0^,^ + a^+v,v)]}- (22) 

In the classical continuum limit we require that: 

Vj + avj -. eg + aaj .„„s 

The interaction term with couplings to external vector and axial vector sources is given by: 
S^ = a^K2 J2 J2('^^^^^,^^^,'"^'-^+'" ~ ^x+vjGl^v^lvi'x)- (24) 

X V 

Our final expression for the lattice action is: 

S = Sg + S^ + S,, (25) 



with Sg, S^ and Sj given in (0), (0) and (|T9|). 

We set Ki = K2 = 1/12 and in the limit of vanishing lattice spacing a -^ identify the 
gauge variables with continuum transport factors through Gx,v = Q^p[igoAfj.{x)v^]. With this 
choice, the action (|25|) has a correct classical continuum limit. 

It is worthwhile noting that our action possesses chiral flavor gauge symmetry under 
transformations: 

Vx -^ [La; h Rx )Vx, 

i^x - ^.(4^ + ^i^^), 
Y -^ (T l^Ji^filU^w r/-t ^-^5 „t 1±75. 

^X,V V 2; r) ~l -'■^X r) ) ■'■ X,V\-'^X+V Q ~l -'■^X+V r\ /' 

Sx + iPx -^ Rx{sx + ipx)Ll, (26) 

which reduces to the form (|^) in the classical continuum limit. 

Generating functional on the lattice can be written as a path integral over the quark and 
gluon fields: 

Z{s,p,v,a)= f DijDijDGe-^^^'^'^'''P''''''\ (27) 



^It has been argued in Ref. |]10| that the use of the BCH lattice actually worsens the problem of lattice 
ferniions as one tries to recover the continuum limit. This is irrelevant for our purposes, because our moti- 
vation for its use comes solely from the need of obtaining better symmetry properties in the strong coupling 
limit. We would not advocate using the BCH lattice in numerical simulations done in the scaling regime, 
where there is good evidence for the recovery of rotational symmetry with the ordinary HC lattice. 



In this expression, s,p,v,a are external sources defined earlier, the action S is given in (^) 
and the integration is performed over all the lattice variables (for example: DG = Y\x,v dGx,v)- 
In the strong coupling expansion, the gauge term (|T7D is suppressed by l/g^. Therefore, 
in the leading order we can neglect it, and the action reduces to S = S^ + Sj. Now the 
gauge variables appear only in S^, moreover they appear only linearly. This implies that the 
product of integrals over gauge variables DG factors into a product of one-link integrals: 

Z = f DipDiPe-^^ n / dG^^^e-^^'^ . (28) 

•^ x,v •' 

Here we have defined: 

Wx,v = -r^{^xi'Gx,vYx,v^x+v - ^|Jx+v^Gl^^Yl^^px). (29) 

One link integrals of this type have been calculated for the case when the gauge field 
couples to a bosonic source in |Q. In ^ it was shown that this result applies equally to an 



interaction of gauge with fermionic fields. The interested reader should consult reference [Q 
for the details. It is important to mention that the bosonization procedure can be extended 
to the next order in l/g^ expansion ^T| and to arbitrary (odd) A^ [|^. Here we just give the 



final result for the generating functional to the leading order in the strong coupling, large-N 
limit: 

Z = j DijD^e''^^,^^'^^^'^^'^^^^-'^ . (30) 

A(x, f) is a matrix in Dirac-flavor space, defined by its element^: 

6 iV 4n _ 

k^.vu = i7^,y: e [^t{^)i^i{^)mi),d]{x+v)iJ]{x+v){i^Yx,.)i,,]. (si) 

F{X) is the following matrix function^ of A: 



l + Vl^A 



F(A) = 1 - v^r^A + In ^ , (32) 

Next we proceed to convert the fermionic path integral into one over bosonic variables as 
in references |Q, ^. In the expression (|30|), A^Tr[F(A)] is the interaction term. It is important 
to notice that it depends only on color singlet combinations ■j/'f (a;)-?/'f (x) at different points 
on the lattice. Since these appear also in the source term 5*^, we can write: 

^3 _ If)- 1 r^S"- 

iVY Nadjf^ "" "^ "^^ Nadjf ^ ' 

Therefore we can rewrite the full generating functional as: 

^Greek letters will denote indices in the direct product of Dirac and flavor spaces, and Latin letters the 
color indices. In this section we will take the number of flavors to be n and only at the end specialize to 
n = 3. 

^F(A) should be understood as a power series in Grassmann variables. 
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with: 

2 4n ^ (9 

The path integral over fermionic variables factorizes into a product of one-site integrals: 

Zo = j Di^Di^e-'^ = n / #,#.e-'^'('^^^'^'^^) = \{ ZoU.), (36) 



zoU.) = j dtp^d^p.e-^ ^^^^^^^\ (37) 

In 1^, p it was shown that Zq{Jx) can be written as an integral over bosonic variables: 

zoU.) = J t/;W,e^^("^'-^^)-^^'^(i'^-^-)+^°"^*. (38) 

Here Ai^ is a unitary bosonic matrix in Dirac-flavor space, the trace is over Dirac-flavor 
indices and the constant is irrelevant so it will be omitted hereafter. 

Using (|3^) and (^5]) we obtain an expression for the full generating functional: 

Z= /■/)>(e^^^-."^'t'^^^=^'"^l+^-^^''^'"-^=^^"^-^''^^°-^=^^^ = f DMe-'^^-^\ (39) 

where A is now a bosonic matrix: 

A(x, v) = ^JY,^,M:,^4YI,M^. (40) 

3o 

The transformation of M. under chiral gauge symmetry follows from (pG]): 

M. -. (L.l^ + fl,i + li)A<,(i.i±25 + flti^). (41) 



One can see that this, together with a law of transformation of Y^^^ (|2^) makes the action 
S{A4) chirally gauge invariant. Now we look for a translationally invariant saddle point of 
the action S{A4) of the form: 

Ml = Mo, (42) 

Uq being proportional to the unit matrix in Dirac-flavor space. From (|40|), using the fact that 
^^ = 1, we obtain Aq = A(x, uq) = Uq/36. For a flavor diagonal mass matrix j^ = Sx = fn and 
no other external sources, the action at this point reduces to: 

7/2 1 1 

S' = -Nj:nFi^) - Y^ln^o + j^amuo]. (43) 



The effective potential is given by: 



By requiring: 



K =-^i^. (44) 

volume 



^ = 0, (45) 

auo 
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we find the saddle point as a positive solution of the quadratic equation: 



«o = ,. , -., (46) 



-22(am) ± ^ {2Aarh)'^ + 1472 
16 + {affi) 

For the massless case Mq = -\/23/2 ~ 2.40. 

We will parameterize the field A^a; in a nonlinear way P: 

M^ = uo exp [^{S^ + 7r,75 + V:^, + zA^7,75 + ^T^a,,)]. (47) 

Here S, vr, V^, A'^ and T'^'' are sixteen Hermitian matrices in flavor space. Of these fields, only 
71 are real Goldstone bosons, as we will show in appendix A. 

In principle, we need to integrate out the scalar, vector, axial vector and tensor fields 
in order to obtain an effective action for pseudoscalar mesons. In the first approximation 
we will neglect these contributions^, and concentrate only on the direct interactions among 
Goldstone bosons:0 

M,, = Uoexp{t^). (48) 

-TO 

Next we expand the action S'(A^) up to fourth order in Taylor series around the vacuum 
( ^2]) and (|i6|). This is sufficient to extract coefficients of effective chiral Lagrangian to order 
p^ since the deviation from the vacuum A^;^^ — Aq is of 0{p). 

The action in (|39| ) can be written as: 

S = -iV5]Tr[F(A,.,,)] - iVnoETr(aj.-M.). (49) 

X,V X 

We expand Tr[F(A2.^„)] to the fourth order: 

Tr[F(A.,.)] = E{- 1^ Tr[(A,,. - Ao)"]}. (50) 

Tl — 1 * \ / Aq 

We evaluate the Dirac traces in a basis where 75 = diag(l, 1,-1,-1). By using the relations 

^!Y{v,a)x^y = Y{v,-a)^^yij, (51) 

which follow from definitions of M. and Y{v,a), we can take the trace over spin indices to 
obtain: 

Tr[(A,,, - Ao)"] = 2A^Tr;[(Z,,,?7i+„iy],,f/, - 1)" + h.c. ]. (52) 

^They are suppressed by the masses of relevant resonances. 

^"in this paper we will not consider the effects of the U(l) anomaly in any detail, so we will restrict our 
attention to the traceless part of tTx ■ 
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Here we have defined new variables which are matrices in flavor space only: 

U^ = exp(z-^), 
-To 

^x,v = 6Xp I [Vx,v + "^x+v.v + {(^x,v + (^x+v,v)\fj 

— tCL 
^x,v ^= exp I - [Vx^v ~r Vx+v,v \(^x,v ~r (lx+v,v)\s- l'-*'^/ 

Tvf indicates a trace over flavor indices, and from now on we will omit the subscript /. 
Evaluating the second term in (|^) is straightforward and one obtains the expression for the 
action in terms of the new variables: 

4 \n /f)'n.fP\ 

S= - 2iV;^^-i -- TT[{Zx,MU.K.Ux-l)-+h.c.] 

x,v n=l '''■ \ "^^ / Ao 

- 2NuoaJ2^Aisx - iPx)Ux + h.c.]. (54) 

X 

Relations (^) , (^) and (|53D determine the effect of chiral flavor symmetry transformations 
on the new variable^: 

Ux — >■ RxUxLx, 

'^ x,v ^ ^x'' x,v-'^x+vi 

^x,v *■ ^x^x,vJ^x+v V-'^j 

The next step is to re-express Tr[(Aa;^„ — Aq)"] in terms of a more convenient variable which 
will facilitate expansion in powers of momentum (or equivalently powers of lattice spacing 
a). The new variable is a, deflned as: 

a = {AMxVi^vUx), 
AMx = Wl^Ux - Ux+.Zl,. (56) 

One can express the traces Tt[{Xx^v — -^o)"] in terms of a: 

Tt[\x,v - Ao] = -2AoTr(a), 

Tr[(A.,, - Ao)2] = 2A2[Tr(a2)-2Tr(a)], 

Tr[(A,,,-Ao)=^] = 2Xl[3TTia^)-TTia% 

Tr[(A,,, - Ao)1 = 2A^[Tr(a^)-4Tr(a3) + 2Tr(«2)]. (57) 

This is very convenient when one truncates the momentum expansion as we will see in a 
moment. We will denote by A^f/^ the coefficients in the expansion of Ayllx in powers of a. 
Relations (|53D and (^) imply that the lowest order term in this expansion is the linear term: 

AMx ^ a{AlUx) + ^{AlUx) + ^{AlUx) + ... (58) 

2 o 



^""^Note that the usual transformation law for U fields follows fi-om the initial form of transformation of 
quark fields (pq ) and the sequence of changes of variables that we performed. 
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Since the expansion of a starts with terms which are quadratic in a, we can neglect terms 
proportional to Tr(«^), Tr(a'') in equations (|57D, and keep only those proportional to Tr(Q;) 
and Tr(a^). 

Next we define the right and left handed external sources analogously to the continuum 
case: 

Any variable which under the lattice symmetry transforms in the same way as the link w"- 
on which it is defined will from now on be denoted by a subscript v. For example: 

ay = -^=-^, a„ = (a„) , etc. (60) 

Partial derivatives in the direction of the link vfj are defined similarly: 

(9^,a = — _ "^ , d^a = dv{dva), etc. (61) 

These derivatives appear in the Taylor series expansion of r^+„ ,; and Rx+v in powers of a: 

rx+v,v = r^ + a(9„r^) + — ((9^r^) + ..., 

Rx+y = R + a{d,R) + —{dlR) + ..., (62) 

and similarly for lx+v,v and L^^+^Q. By substituting the previous Taylor series expansions 
into equations (|53D and (^) one can derive the laws of transformation of r^ and /^: 

r„ -^ Rr^R^ + iR{dyR^), 

/„ -^ Ll^L^ +iL{dM)- (63) 

Since this is true for any vf^ (pOD, in the classical continuum limit (p3|) this law of transfor- 
mation becomes the usual law of transformation of gauge fields (|^). Let us now define D^U 
on the link v in analogy with the continuum covariant derivative: 

D^U = d.,U-ir^U + iUl^. (64) 

In the classical continuum limit, D^U reduces to its continuum counterpart (p!2|). More- 



over, from the transformation laws (33), one can see that D^U transforms covariantly on the 
lattice: 

D„U -^ R{DM)LK (65) 

We now expand A^Ux in a Taylor series in a by a straightforward but tedious calculation and 
then re-express coefficients in this expansion, defined in (^8]), in terms of the lattice covariant 
derivatives Z)„^:0 



Apx = -Dp-2tr,DM, 



Apx = - Dp ~3tr,DtU- Slider,) D,U + SriD,U. (66) 



"'^■^Here we are using the convention that, when there is no subscript x, a variable (or its derivative), is 
evaluated at x. 

^^Here by DlU we mean D„{D„U) etc. 
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Although not all coefficients in the expansion of AyU are covariant under chiral transforma- 
tions, Tr(Q;) and Tr(a^) are manifestly covariant! 

Tr(«) = TT{a\DMy{DM) + a^[ip.f/)T(Z},3f/) + h.c. + ]{Dluy{DlU)]}, (67) 

o 4 

Tr(a2) =Tr[a\DMy{DM){D^Uy{DM)]- (68) 

In appendix B, we derive results for the summation of tensors of rank two and four over 
the twelve links in the positive direction: 

^aA = 3^a^6^ (69) 

Y^ a^KcJ^ = - Y,{a^h^'c^(r + ai.c^b^d" + a^d^^Wc^"). (70) 

By using these expressions, summing over links and Wick rotating back to Minkowski space, 
we get manifestly Lorentz and gauge invariant expressions for the traces: 

^Tr(a) = Ti{^[{D^U)\D^D,D''U) + {D^U)\D''D^D,U) + {D^Uy{D-D,D^U) 

o 

ETr(a') = lTT[{D^UyiD^U)iDM)HD''U) + iD^U)\DM)iD^'Uy{D''U) 

^{DJiJ)\D,U){D-U)\D^U)\. (71) 

The last step is to substitute these expressions back into (|58D . The action on the lattice gives 
the Lagrangian after dividing with the volume of space-time on the BCH lattice: 



a^ 



E = ^ld'x (72) 



We can transform all the terms of order p*^ into the canonical basis of ^ using the property 
of SU(3) matrices derived there: 

Ti[iD^UyiDM)iD^UyiD''U)] = -2Tr[(D^f/)t(D^f/)(D'^f/)t(D,f/)] 

+ i{Tr[(D^'t/)t(D^f/)]}2 + TT[{D^uy{D''U)]TT[{D^Uy{DM)], (73) 

as well as the equation of motion at order p^: 

{D^'D^U)U^ - U{D>^D,Uy = xU^ - Ux^ - ^Tt{xU^ - Ux^). (74) 

We have collected these results in Table 1. and compared them with the coefficients which 
have been determined experimentally at the scale of rrifj in Q] 
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theory 


experiment 


Li 


0.5 (Foa)2 


0.9 


± 0.3 


L2 


1.1 (Foay 


1.7 


± 0.7 


L3 


10.4 {Foaf 


-4.4 


± 2.5 


U 


(Foa)2 





± 0.5 


L5 


(Foa)2 


2.2 


± 0.5 


u 


(Foa)2 





± 0.3 


Lr 


0.9 (Foa)2 


-0.4 


± 0.15 


Ls 


2.6 (Foa)2 


1.1 


± 0.3 


L9 


-6.7 (Foa)2 


7.4 


± 0.2 


Lio 


-6.7 (Foa)2 


-5.7 


± 0.3 



Table 1: Comparison of the strong coupling expansion with experimental results from 01, [0]- 
The left column contains our results from the strong coupling expansion. The right column 



contains the experimental values L[ 



rrir 



All numbers are in units of 10" 



4 Beyond the strong coupling regime 

We have chosen to express the entries in Table 1. in terms of a common factor (-Foa)^ 
because of a point we will be making below, but the value of this quantity is also determined 
by the strong coupling expansion, from the 0{p^) terms in the lattice chiral Lagrangian, 
and is given by 2.1A^. With 3 colors the entries in the table come out about one order of 
magnitude higher than the values derived from experiment. We take this to be an indication 
that the strong coupling limit produces too tight a binding between quark and antiquark, 
which in turn leads to an unacceptably high value when Fq is expressed in terms of a~^ 
(or, equivalently, to too large a value for a). If we allowed ourselves to set a~^ ~ Fq, 
which would be a more reasonable scale for a strong coupling calculation, the theoretical 
predictions for several coefficients would compare rather well with the experimental results. 
Of course we should not make too much of this agreement, because we expect the values of 
the coefficients to depend on detailed dynamical features of the interactions which may not 
be well reproduced by the strong coupling approximation. What is much more important 
is the demonstration, through the strong coupling expansion, of the main point we wanted 
to make, namely that one can formulate an effective chiral Lagrangian on the lattice and 
that this can be the vehicle for the derivation of the parameters of the continuum chiral 
Lagrangian. 

Another point that should be mentioned is that the expansion for small lattice momenta 
that we have used to derive the coefficients corresponds to using the lattice chiral Lagrangian 
in the tree approximation. We have calculated a few of the corrections that would be induced 
by one loop diagrams and have found these to be small. However, the same factor of (Foa)"^ 
which we discussed above enters also as a coefficient in the denominator of the loop diagrams, 
so that the statement that the loop corrections are small should be taken with caution. 
They would become larger if one could assign to Fqu the smaller value which we advocated 
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above. Nevertheless, even if the loop corrections turned out to be large in a more realistic 
calculation, this would only represent a technical problem and not a conceptual difficulty 
for the whole approach. The very important fact is that, because we are working with 
a regularized theory, all physical quantities are well defined and finite. Moreover, since the 
lattice system also exhibits spontaneous breaking of chiral symmetry and Goldstone bosons in 
the limit of vanishing quark mass, it will always be possible to perform an expansion for small 
lattice momenta. The coefficients of such expansion are in any case well defined quantities. 
Thus the question is whether they can be calculated by a perturbative expansion of the chiral 
effective theory, which would be more convenient, or whether, failing such possibility, one 
will have to resort to numerical techniques. Nevertheless, even in the latter case, one would 
still be dealing with a bosonic system, and would therefore avoid the need of simulating the 
quantum fluctuations of fermionic variables, which represents today the major difficulty one 
faces in the implementation of QCD numerical simulations. 

Finally, in order to obtain reliable values for the coefficients one should perform the 
calculation in the range of values of the lattice coupling constant (/3 ~ 5.7 — 6, including 
quark degrees of freedom in the simulation) where one witnesses the onset of scaling to the 
continuum. Such a calculation can at present only be done by numerical techniques. The way 
we envisage it could be carried out would be by assuming a sufficiently large set of couplings 
for the lattice effective theory and fixing them through the matching of an overcomplete 
set of expectation values. This is very similar to procedures commonly used in numerical 
studies of renormalization group transformations. The crucial point is that, because the 
lattice effective theory already accounts for the long range excitations of QCD, the matching 
should only require a reasonably small lattice size, of the order of the inverse p mass. The 
two theories (the one derived from the original QCD Lagrangian and the effective theory) 
should produce exactly the same values for the observables on any lattice size, because they 
are mathematically equivalent. Working with a small lattice, we are confident, one would be 
able to fix the parameters of the lattice effective theory with a good degree of accuracy and 
would thus establish a solid base for the subsequent determination of the parameters of the 
continuum effective chiral Lagrangian. 
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Appendix A. Goldstone bosons 

In this appendix, we will show that the fields tt^ are Goldstone bosons, corresponding to 
the spontaneous breaking of the chiral symmetry. 

It was shown in P, 0, for the case of no external sources on the hypercubical lattice, 
that the action in (|39|) exhibits spontaneous chiral symmetry breaking: 

U{n) X U{n) -^ U{n). (75) 

In that context, n represented the number of color triplets of staggered fermions. Our situ- 
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ation is somewhat different: we start from full-ffedged fermions with Dirac, ffavor and color 
indices and n represents the number of flavors. Moreover, we have a different lattice structure. 
We obtain the inverse propagator in moment umf^ space for n^ fields tTj, from expressions 
and ( |57| ) by expanding to the quadratic order in powers of fields tt^ : 

D{p) = B[l2-Y,cos{p-v)], (76) 

V 

AN 
5 = -t[AoF'(Ao) + A^F"(Ao)]. (77) 



Fr 



We can rewrite the momentum dependent part of D{p) in terms of momentum components 

> cos(»-f)= 2 cos ^^(cos ^^ + cos ^^ + cos ^^) + 

V ^ x/2^ 72 y/2 ^' 

COS ^^(cos ^^ + COS ^^) + COS ^^ COS ^^ . (78) 

v^^ V2 72 72 72^ ^ ^ 



By substituting this back into the expression ( [76D one can immediately see that the propagator 
D~^{p) has a massless pole p^ = im = at the origin p = for each of the n^ fields contained 



in the matrix iix- 



Appendix B. Lorentz covariance 

In this appendix we derive relation (p5|). From the definition (pDj) 



^ _ M-^ ai + attj hi + ahj 

2_^ ttyOv — 2^ -J= -J= 

= 2 ^ ^'^*^* "^ ^^^^ ^ a(ai&j + aj6i)]. (79) 

When summing over a = ±1, terms proportional to a cancel and the others get a factor of 
2: 

^ a^hy = ^{aibi + ajbj). (80) 

Finally, we perform the sum over I < i < j < A and obtain: 

^aA=3^a^6^ (81) 

Relation ( [70D can be proven in the same way. 
^■^In this appendix, for simplicity we take a = 1. 
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